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Abstract — This work considers optimizing LDPC codes in 
the physical-layer network coded two-way relay channel using 
noncoherent FSK modulation. The error-rate performance of 
channel decoding at the relay node during the multiple-access 
phase was improved through EXIT-hased optimization of Tanner 
graph variable node degree distributions. Codes drawn from 
the DVB-S2 and WiMAX standards were used as a basis 
for design and performance comparison. The computational 
complexity characteristics of the standard codes were preserved 
in the optimized codes by maintaining the extended irregular 
repeat-accumulate (eIRA). The relay receiver performance was 
optimized considering two modulation orders M = {4,8} 
using iterative decoding in which the decoder and demodulator 
refine channel estimates by exchanging information. The code 
optimization procedure yielded unique optimized codes for each 
case of modulation order and available channel state information. 
Performance of the standard and optimized codes were measured 
using Monte Carlo simulation in the flat Rayleigh fading channel, 
and error rate improvements up to 1.2 dB are demonstrated 
depending on system parameters. 

I. Introduction 

The two-way relay channel (TWRC) models two source 
nodes outside radio range exchanging information through a 
single relay node in range of both. Physical-layer network 
coding (PNC) |[T| may applied in the TWRC to increase 
throughput by reducing the number of time slots required for 
exchange. The exchange between the sources is broken into 
two phases; the multiple-access (MA) phase and broadcast 
phase. During MA, the sources transmit at the same time 
to the relay, and the signal received at the relay is the 
electromagnetic sum of signals transmitted by the sources. 
During the broadcast phase, the relay transmits the network- 
coded sum of signals to the tranceivers, and each performs 
network decoding to recover the information transmitted by 
the opposite sources. 

Previous work developed a demodulation and decoding 
scheme for the relay in the PNC-coded TWRC using M-ary 
FSK modulation Q ||^. This scheme supports iterative decod¬ 
ing between decoder and demodulator to improve error-rate 
performance of the received network-coded information bits. 
The receiver supports power-of-two modulation orders and 
improved decoding performance based on available channel 
state information (CSI). A simulation regime was performed 
to demonstrate the performance of the scheme using Turbo 
decoding at the relay. 

The error-rate performance of LDPC codes depends on the 
structure of the parity check matrix and thus the properties 
of the Tanner Graph and degree distribution. In this work 
we apply an optimization technique which identifies variable 
node degree distributions likely to yield good codes using ex¬ 
trinsic information tranfer characteristic (EXIT) curve fits Q . 
LDPC codes taken from the DVB-S2 and WiMAX standards 
are used as a basis for design and performance comparison, 
and their variable node degree distributions are optimized to 
yield codes exhibiting error-rate performance superior to the 


standard codes. The complexity characteristics of the standard 
codes are maintained in the optimized codes by preserving the 
extended irregular repeat-accumulate constraint. 

A variety of optimization techniques exist for improving 
LDPC code performance. Random permutation matrices, large 
girth optimization, and column weight optimization have been 
applied to WiMax standardized codes yielding gains up to 
0.4 dB Q. Joint design of and LDPC codes under carrier 
frequency offsets demonstrate the relationship between system 
properties and code performance |[^. 

A PNC scheme compensating for correlation among trans¬ 
mitted symbols introduced by channel coding to improve 
decoder performance at the relay has been considered Q. 
In this work it is demonstrated that LDPC codes having 
different degree distributions are preferred in different SNR 
regimes. Roughly, smaller degrees are preferable at low SNR 
and higher degrees at high SNR. A PNC system used as a basis 
for optimizing LDPC codes via EXIT analysis demonstrated 
performance 1.5 dB from capacity |j^. 

The rest of this work is organize as follows. Section II 
develops the system model applied for code optimization and 
error-rate simulation. Section III describes the variable node 
degree distribution optimization technique in information- 
theoretic terms, while Section IV describes the computational 
procedure followed to generate the optimized codes and 
demonstrates the performance of the codes versus standard. 
Section V provides concluding remarks. 


IT System Model 

This section describes the model applied for LDPC code 
optimization and error-rate simulation. The model is depicted 
in Eigure [T| 


A. Transmission by Source Nodes 

The source nodes Afi,i G {1,2} generate binary infor¬ 
mation sequences = [ui^i, ...,UK,i] having length K. A 
rate-i? LDPC code is applied to each u^, generating a length 
N — K/R channel codeword, denoted by b' = [5i i..., 6^ 
The codeword is passed through an interleaver, modeled as a 
permutation matrix 11 having dimensionality N x N : hi = 
b'n. Let V = (0,..., M — 1} denote the set of integer indices 
corresponding to each ESK tone, where M is the modulation 
order. The number of bits per symbol is p, = log, M. The 
codeword b^ at each node is divided into Ng = N/p sets of 
bits, each of which is mapped to an M-ary symbol G V, 
where k denotes the symbol number, and i denotes the source. 

The modulated signal transmitted by source A/} during 
signaling interval kTg < t < {k -\- l)Ts is 





2tt 



{t - kTs) 


( 1 ) 


where is the transmitted signal, A is the carrier 

frequency, and Ts is the symbol period. 

The continuous-time signals Sk,i{t) are represented in 
discrete-time by the set of column vectors ^ : qk,i G T)}. 
The column vector . is length M, contains a 1 at vector 





Fig. 1. System Model - TWRC DNC MAC Phase 


position ( 7 fc i, and 0 elsewhere. The modulated codeword 
from source Ni is represented by the matrix of symbols 
Xi = i], having dimensionality M x Ng, where 

.. 


B. Channel Model 

All channels are modeled as flat-fading channels having 
independent gains for every symbol period. The gain from 
node Ni to the relay during a particular signaling interval k is 
denoted by hk^i- The gain is represented as hk,i = a:k,ie^^'°'\ 
where ^ is the received amplitude and 6k,i is the phase. The 
received amplitude has a Rayleigh distribution with parameter 


(T = y and 6k,i is the phase having a Uniform distribution 
over interval [0, 27r). 

Consider transmission of a single frame of Ng symbols to 
the relay. The received frame is 


Y = XiHi -b X2H2 + N ( 2 ) 

where is an Ng x Ng diagonal matrix of channel co¬ 
efficients having value hk,i at matrix entry {n,n),n G 
{1,2,Ng} and 0 elsewhere, and N is an M x fVq noise 
matrix. A single column of Y represents a single signaling 
interval, is denoted by y, and referred to as a channel obser¬ 
vation. In terms of this definition, Y = [yi,..., yAr ], where y^ 
denotes the /c-th channel observation. Denote the ^-th column 
of N by rifc. Each column is composed of zero-mean circu¬ 
larly symmetric complex Gaussian random variables having 
covariance matrix NqIm\ i-S-, ~ Nc{0, NqIm)- Nq is 

the one-sided noise spectral density, and Im is the M-by- 
M identity matrix. The noise spectral density is defined as 
-^0 = lox/iouM ’ where X is the bit-energy to noise ratio in 
decibels Increasing modulation order decreases the noise 
spectral density to reflect an increase in energy per-symbol 
when utilizing constant energy per-bit. 

C. Relay Reception 

The goal of the digital-network network coding (DNC) 
relay receiver is to detect the network-coded combination of 
information bits transmitted by the end nodes, u = Ui 0 U 2 . 
The relay receiver takes as input the frame of channel observa¬ 
tions Y. The symbols transmitted by the sources are assumed 
to be received in perfect synchronization at the receiver. 
Demodulation and iterative channel-decoding are applied to 
detect u. Define the network codeword as 


b = bi © b2 

= [ ^o(x/c,i) © bo{Xk,2) ■■■ (>^-l(Xfc,i) © 6^_i(Xfe_2) ] 

V fee iVq} (3) 

where bm{^k,i) denotes the m-th bit mapped to the fc-th 
symbol transmitted by source N- Since the channel code 
considered is a systematic linear code, b forms a code from 
the codebooks used by the source nodes, thus, the channel 
decoding operation yields a hard decision on the network 
coded message bits u. 

The demodulator is implemented using super-symbol prob¬ 
ability mapping and DNC soft bit-mapping (SOMAP), having 
formulation described in ©■ 

The super-symbol probability mapper takes as input the 
matrix of received symbols Y and produces estimates of the 
probability of receiving each super symbol qk, defined as the 
tuple 


Qk = {qk,i,qk, 2 ) qk,i,qk, 2 &'D qkGDxV (4) 

to produce estimates P{qk), 0 < k < Ng. The cardinality 
of D X D is M^, thus, the number of probability estimates 
computed at the relay during each symbol period grows 
exponentially with modulation order. This computational com¬ 
plexity is a practical consideration for implementing the relay 
receiver. 

The super-symbol probability estimates are passed to the 
DNC SOMAP along with with extrinsic decoder feedback v 
to produce extrinsic information z, likelihood ratios of the 
network-coded bits b communicated by the frame. 

The extrinsic information is deinterleaved to produce z' = 
zll”^ and passed to the decoder. The decoder refines the es¬ 
timate of z', producing new extrinsic information v' which is 
interleaved to produce v = v'll and returned to the SOMAP. 
After a specified number of iterations, a hard decision is made 
on the network coded bits u. 

The transmitted symbol energies and noise spectral density 
Nq are known at the demodulator. The demodulator may 
operate under several cases of channel state information (CSI) 
©; the case in which the gain are completely known (full 
CSI), the case in which only the fading amplitudes ak,i are 
known (partial CSI), and the case in which no information 
about the gains is known (no CSI). 

Since the focus of this work is on the relay reception phase, 
details of the relay-to-source broadcast phase are omitted. 
Interested readers are referred to ©. 

III. LDPC Code Optimization 

This section describes the application of a procedure for 
optimizing LDPC code error-rate performance based on © 
under a range of receiver configurations to the DNC relay 
reception phase by varying the LDPC code degree distribu¬ 
tion. Relevant LDPC concepts are introduced, followed by a 
description of metrics relevant to optimization. A theoretical 
description of the optimization technique is provided followed 
by the procedure t^en. 

The goal is to produce LDPC code designs having the 
same rate and complexity as off-the-shelf codes from the 
DVB-S2 and WiMAX standards, while improving error rate 
performance through optimization of the parity check matrix. 
The codes are constrained to the class of extended irregular 
repeat accumulate (eIRA) codes. This constraint ensures com¬ 
putational efficiency by guaranteeing systematic encoding © . 
The codes considered are check-regular, meaning that every 
check node has the same degree, dc. 

















































An LDPC code is a linear block code having a sparse parity 
check matrix H, which is an N — K x N binary matrix. 
An alternative representation for the parity check matrix is 
the Tanner Graph, a bipartite graph in which one of the sets 
contains variable nodes, and the other contains check nodes 
pO) . Each column of H corresponds to a variable node, and 
each row a check node ng- A 1 entry in H at row and column 
{k, n) corresponds to an edge between the fc-th check node 
and n-th variable node. The degree of a node is its number 
of incident edges, thus, the number of I’s in each row and 
column for each check and variable node, respectively. 

The variable and check nodes may be modeled as soft-input, 
soft-output decoders denoted as VND and CND, respectively 
0. The VND decoder takes a channel LLR L^yND and a- 
priori information La,vND as input and produces a-posteriori 
information Lo,vnd as output. The extrinsic information is 
defined as LeyND = LayND — LoyND- The mutual infor¬ 
mation between the L^-yMo and the corresponding network- 
coded bit is IsyND, while the mutual information between 
LayND and the network-coded bit is lAyND- 

The CND decoder takes a-priori information La,CND as 
input and produces a-posteriori information Lo,cnd as output. 
The extrinsic information is Le,CND = La,CND — Lo,cnd- 
The mutual information between the Le,CND and the cor¬ 
responding network-coded bit is Ie,cnd, while the mutual 
information between La,CND and the network-coded bit is 
Ia,cnd- 

An EXIT chart visualizes the transfer characteristic of a 
decoder by plotting the a-priori mutual information against 
the extrinsic. It has been shown that matching the VND 
and CND transfer characteristics of a particular LDPC code 
through selection of variable node degrees improves error rate 
performance HD- 


A. Optimization through Selection of Variable Node Degree 

This section describes the analytic formulation for degree 
optimization. The degree distribution is defined as the set 
containing the degrees assigned to the variable nodes, and the 
number of variable nodes taking each particular degree. De¬ 
note a degree distribution as V = { : oi,..., dy,D '■ od }, 

where Oi is the number of variable nodes having degree dy^i. 
For a particular code, all check node degrees are the same, 
thus, the code is check-regular. Denote the check node degree 
for a particular code as dc- Define as the fraction of variable 
nodes having degree dy^i 


a, = ^, i = (5) 

each taking a value 0 < Oi < 1, and = 1- 

A realizable degree distribution satisfies the constraints 
imposed by the LDPC code parameters. The sum total of edges 
incident on the variable and check nodes must be equal. The 
number of edges incident on variable nodes having degree i 
is Cy^i = Oidy^iN, thus, the total number of edges incident on 
all variable nodes is 

D D 

Cy - ^ ^ ^ ^ O^idy^iN ( 6 ) 

and the total number ofiedges incident on the check nodes is 

ee = daiN - K). (7) 


Equating e„ and Cc and rearranging, 

Oidy^iN _ ^ 


( 8 ) 


^^daiN-K) 

A particular degree distribution V is selected by choosing 
parameters V and dc which satisfy and the sum constraint 


on Qi- The challenge is to select a V which optimizes error 
rate performance for particular channel conditions and relay 
receiver configurations. 

The optimization procedure requires generating EXIT 
curves corresponding to the variable-node decoder (VND) and 
check-node decoder (CND). The VND curve is described by 
the variable node degree distribution V and the mutual infor¬ 
mation computed between the DNC-SOMAP output and the 
network coded bits. The CND curve is completely described 
by the degree of the check nodes dc, and is fixed for all 
variable node degree distributions. 

Generation of the VND curve begins by computing the 
mutual information between the soft-values at the output of 
the DNC-SOMAP and the network-coded bits, denoted as the 
detector characteristic Ie,det(.Ia,deT: ^)- Ia.det is the 
mutual information between the DNC-SOlviAP a-priori input 
and the network-coded bits, and ^ is the signal-to-noise ratio 
of the super-symbol. The detector characteristic Ie,det is 
generated by Monte Carlo simulation under the assumption 
that the a-priori input Ia.det is conditionally Gaussian. 
The Gaussian assumption can be verified empirically using 
histograms of several decoding runs HI- 

The VND curve is found by computing a separate curve 
for each variable node degree and combining each to form 
the final curve describing the entire variable node detector. 
The curve corresponding to the i-th degree is given by 


lEyND{lAyND,dy^i,lE,DET) — ■■■ 

J (^\J(dy.i — l)J~^{lAyND)^ + [J~^{Ie.detW'^ (9) 


where the J-function is defined in |13| and is c omp uted 
using the truncated-series representation presented in ]14[ , and 
lAyND = Ia,det- The VND curve describing the entire 
detector is given by 


lEyND{lAyND, Ie.det) = ■■■ 

D 

bi ■ lEyND{lAyND,dy^i, Ie.det) (10) 

i=l 

where bi — is the fraction of edges incident 

on variable nodes of degree i. 

The CND curve is computed according to Q 

Ie,cnd{Ia,cnd, dc) = I — J (^\/dc — I ■ J~^{1 — Ia.cnd) 

( 11 ) 

where Ia.cnd is the mutual information between the check 
node inputs and the network-coded bits. Note that the CND 
curve is independent of the particular variable node degree 
distribution V. 

The curve fit for the VND and CND mutual information 
metrics corresponding to a particular degree distribution is 
performed by plotting the curves for a range of values of ^ 
and noting the value in which the curves barely touch Q. 
Specifically, the VND curve is plotted with lAyND on ffie 
horizontal axis and Ie.vnd along the vertical. The CND curve 
is plotted with Ie.cnd along the horizontal and Ia.cnd 
along the vertical. The value of ^ is varied from highest 
specified point to lowest, and the liighest point for which the 
curves touch is defined as the EXIT threshold. Considering 
a particular relay receiver and channel configuration, the 
optimal variable node degree distribution is considered as the 
distribution having the lowest EXIT threshold. 






IV. EXIT-Optimized LDPC Code Pereormance 

This section presents the error-rate performance of EXIT- 
optimized LDPC codes and the simulation optimization proce¬ 
dure used to generate the codes. Several cases of modulation 
order and relay receiver CSI are considered. LDPC codes 
defined by the DVB-S2 and WiMax standards are used as 
reference codes. 

A. Optimization Procedure 

The purpose of this subsection is to describe the procedure 
for generating EXIT curves describing the performance of 
selected degree distributions. 

The procedure involves determining the detector EXIT 
characteristic through simulation for a chosen receiver con¬ 
figuration and channel model, and using this characteristic to 
generate the combined variable node/detector and check node 
characteristics over a range of SNR Q. 

The notation [•] following a bold variable denotes a vec¬ 
tor element, for instance, Z[fc] denotes the fc-th element of 
one-dimensional vector Z and W[i, k] denotes the (i, fc)-th 
element of two-dimensional vector W. 

1) Detector Characteristic: This section describes the com¬ 
putation of the detector characteristic for the DNC relay re¬ 
ceiver. These values are computed via Monte Carlo simulation 
as follows 

1) Select source and relay frame size L, modulation order 
M, received channel state information (CSI or no CSI) 
and SNR value. 

2) Select a discrete range of a-priori demodulator mu¬ 
tual information values, represented as the vector 
^A.DETlk] — k/B,0 < k < B. In all of our exper¬ 
iments, the value B = 100 provided visibly smooth 
detector transfer characteristics. 

3) Simulate a frame of channel-corrupted relay-received 
symbols Y according to the model given in Section 
having length Nq = L/p and modulation order M. 

4) Eor every value of 

• Generate a vector of length L a-priori LLRs v 
representing LDPC decoder feedback under the 
assumption that the feedback is conditionally Gaus¬ 
sian 

V = (b - l/2)crfe -f XCT/c (12) 

where cr^ = J~^(iA,DET[k]) Q and x is a length- 
L vector of i.i.d. samples or i standard normal 
distribution. 

• Compute the decoder extrinsic information z ac¬ 
cording to the receiver formulation provided in 1^ . 

• Compute the simulated value of lE,DET[k] as fnl 

lE^DETik] = 1- -^-y^^max*(0,z[£](-l)‘’M) 
° ^ £=0 

(13) 

• Numerically fit a third-order polynomial to the 
sampled detector characteristic using Ia.det as 
the abscissa and 1e,det as the ordinate, yielding 
coefficients /o ,...,/3 

lE,DET[k] =/3Ia. DET[k]^ + / 2 IA ,DET[kf + 

fliA,DET[k] + /o- (14) 

Define function foET^x) = fsx^ + f 2 x'^ + fix +fo. 
This function approximates the detector characteris¬ 
tic and is used in the computation of the combined 
detector/decoder characteristic. 


2) Combined Detector/Decoder Characteristic: This sub¬ 
section describes the simulation procedure for generating the 
combined detector/decoder characteristics. 


1) Select a discrete range of a-priori demodulator mu¬ 
tual information values, represented as the vector 
'^AyND[k\ = k/B,0 < k < B. In all of our exper¬ 
iments, the value B = 100 provided visibly smooth 
detector transfer characteristics. 

2) Eor every value of detector a-priori mutual information 
iA,VND, 

• Eor every degree dy^i compute the a-priori informa¬ 
tion input to the detector according to 

I-A.DEtI*) k] = J dy^i ■ J • 

(15) 


• Compute lEyND according to (lOi 


D 

'^EyND[k] = Y,h-j[y/G) (16) 


where 


G={dy^i — l)J ^(lA,yAr£)[fc]) + ■■■ 

[J-\fDET{l'A,DET[tM)?- ( 17 ) 


The check node characteristic is completely determined by 
the parameters of the code and is not dependent on the prop¬ 
erties of the detector output. Noting that Ie,cnd = I a.vn d 
compute the a-priori CND characteristic according to (111 

— "iE,CND[k]) 


l-A,CND 


[k] = J 


Vdc - 1 


(18) 


Base 

Standard 

M 

CSI 

y 

^ fit 

JVn 

DVB-S2 

4 

Partial 

El = { 2:25920, 4:34560, 22:4320 } 

11.9 




V 2 ={ 2:25920, 33:4560, 30:4320 } 

12 




V 3 ={ 2:25920, 4:37152, 49:1728 } 

12 



None 

V 4 ={ 2:25920, 3:34560, 30:4320 } 

12.3 




Va, ={ 2:25920, 4:34560, 22:4320 } 

12.3 




Ve ={ 2:25920, 4:37152, 49:1728 } 

12.4 


8 

Partial 

Er ={ 2:25920, 3:34560, 30:4320 } 

9.5 




Eg ={ 2:25920, 35:3650, 3:35230 } 

9.5 




Eg ={ 2:25920, 4:34560, 22:4320 } 

9.7 



None 

Eio ={ 2:25920, 3:34560, 30:4320 } 

9.8 




Ell ={ 2:25920, 35:3650, 3:35230 } 

9.9 




Ei2 ={ 2:25920, 4:34992, 24:3888 } 

10.1 

WiMax 

4 

Partial 

Ei3 ={ 2:672, 3:96, 3:1296, 9:240 } 

12.9 




Ei4 ={ 2:672, 3:96, 3:1356, 11:180 } 

12.9 




Eis ={ 2:672, 3:96, 3:1376, 12:160 } 

12.9 



None 

Eie ={ 2:672, 3:96, 3:1248, 8:288 } 

13.3 




Ei7 ={ 2:672, 3:96, 3:1296, 9:240 } 

13.3 




Eis ={ 2:672, 3:96, 3:1356, 11:180 } 

13.3 


8 

Partial 

Ei9 ={ 2:672, 3:96, 3:1296, 9:240 } 

10.4 




E 20 ={ 2:672, 3:96, 3:1356, 11:180 } 

10.4 




E 21 ={ 2:672, 3:96, 3:1376, 12:160 } 

10.4 



None 

E 22 ={ 2:672, 3:96, 3:1356, 11:180 } 

10.8 




E 23 ={ 2:672, 3:96, 3:1376, 12:160 } 

10.8 




E 24 ={ 2:672, 3:96, 3:1392, 13:144 } 

10.8 


TABLE I 

Degree Optimization Results 


B. Optimization Results 

In this subsection we apply the optimization procedure to 
improve performance of standardized DVB-S2 and WiMax 
LDPC codes. The variable node degree distributions for two 
standardized codes are optimized, and their error rate per¬ 
formance is compared under simulation in the DNC relay 
reception phase. 





















Fig. 2. Example EXIT fit - DVB-S2 constraint 


The first code is taken from the DVB-S2 standard |T^, 
having codeword length N = 64800 and rate R = 3/5. The 
second code is specified by the WiMax standard having 
length N = 2304 and rate R = 2/3. The parity checK matrices 
for these codes are specified in their respective standards. Note 
that the WiMax standard specifies two distinct parity check 
matrices for rate 2/3 codes, denoted as 2/3A and 2/3B. We 
consider case 2/3A. These codes satisfy the eIRA constraint 
and are check-regular, with dc = 11 for the DVB-S2 code and 
dc = 10 for WiMax. 

The parity check matrix for each standard codes con¬ 
tains a sub-matrix satisfying the extended irregular repeat- 
accumulate (eIRA) constraint. This constraint simplifies en¬ 
coding and decoding complexity. The submatrices comprising 
H are 

H=[Hi|H2] (19) 

where H 2 satisfies the eIRA constraint and has (N — K) rows 
and {N — K) columns. To preserve the complexity advantages 
of eIRA, we retain H 2 exactly as specified by the standards 
and consider optimizating the variable node degrees specifying 



Fig. 3. Bit Error Rate - DVB-S2 constrained, EXIT-optimized codes, M=4 


Hi. Retaining H 2 places constraints on V such that codes 
based on DVB-S2 have dv,i = 2, oi = 25920, and codes based 
on WiMax have dy,i = 2,oi = 672 and 2 = 3, 02 = 96. 
All other degrees may be chosen freely. 

Variable node degree optimization is performed under sev¬ 
eral cases of receiver configuration and channel state infor¬ 
mation, specifically, modulation orders M = 4, 8 with and 
without CSI, for a total of four configurations. A range of 
degree distributions is considered for each code and receiver 
configuration, and the best performing degree distribution of 
each is realized and simulated. 

Considering the DVB-S2 LDPC code, the number of dis¬ 
tinct variable node degrees is D — 3, and the degrees 
considered are all combinations of the sets dy^i = 2, 
dv ,2 G {3,4,..., 100}, and dy ^3 G {dv ^2 + Ij dy ^2 + 2,..., 100} 
which satis fy the constraints for realizable codes described in 
Subsection |III-A| Considering the WiMax code, the number 
of distinct variable node degrees is = 4, and the degrees 
considered are all possible combinations of dy i = 2, dy 2 = 3, 
d „,3 G {1, 2,..., 100}, and G {1, 2,..., 100}. 

For every receiver configuration and degree distribution V, 
the EXIT threshold is computed to discover the the single 
variable node degree distribution having the best performance. 
For each distribution, the threshold is defined as the highest 
SNR value ^ for which the VND and CND decoder EXIT 
curves touch. The best performing distribution is defined as 
minimizing this SNR value. An example is shown in Fig. S 

The optimization results are shown in Table |I] For each 
degree distribution V, a code is realized having parity check 
matrix satisfying the degree distribution and maintaining the 
eIRA constraint. The parity check matrix is generated ran¬ 
domly. In the interest of space, further details of parity check 
matrix generation are omitted here. 

C. Error-rate Performance of Optimized Codes 

The error-rate performance of the optimal randomly gener¬ 
ated codes is compared to the standardized codes via Monte 
Carlo simulation according to the model described in Section 

m 

The source nodes generate information sequences having 
length K — 38800 for the DVB-S2 codes and K — 1536 for 
WiMax. The information sequences are channel encoded to 
produce codewords having length N = 64800 for DVB-S2 


M=8 

N=64800 

K=38880 

rate=3/5 

iterations=100 

Solid lines; partial CSI 
Dashed lines: no CSI 


-Standard, BICM 

Standard, BICM-ID 


-Standard, BICM 
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Fig. 4. Bit Error Rate - DVB-S2 constrained, EXIT-optimized codes, M=8 

















































M=4 


N=2304 

K=1536 

rate=2/3 

iterations=100 


Solid lines; partial CSI 
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no CSI 
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Fig. 5. Bit Error Rate - WiMax constrained, EXIT-optimized codes, M=4 



Fig. 6. Bit Error Rate - WiMax constrained, EXIT-optimized codes, M=8 


and N = 2304 for WiMax. These codewords are mapped to 
FSK symbols. Modulation orders M = {4,8} are considered. 
The modulated symbols are transmitted over the channel and 
received at the relay. The relay performs decoding considering 
the cases of having channel state information and no channel 
state information. 

In all cases the relay receiver performs 100 decoding itera¬ 
tions before making a hard decision on the network-coded bits. 
The standard codes are simulated for two cases of decoder-to- 
demodulator extrinsic information feedback: BICM, in which 
no information is fed back from decoder to demodulator, and 
BICM-ID, in which information is fed back to the demodulator 
after each decoding iteration. Simulation of the randomly 
generated codes always uses BICM-ID. 

Error rate performance for DVB-S2 and WiMax codes, 
M =^4, 8T and for partial and no CSI are shown in Figures 
l^ffll^andlm The degree distributions for each plotted code are 
aenoted onme legends as Vi, referring to a degree distribution 
listed in the column V in Table |I] 


Consider the results for DVB-S2. For M = 4 and partial 
CSI, the best performing optimized code shows an energy 

efficiency improvement of approximately 1 dB, while the no 
CSI case shows improvement of roughly 0.8 dB. For M = 8 
and partial CSI, the improvement is about 1.2 dB, while no 
CSI shows 1 dB of improvement. 

Now consider WiMax results. For M = 4 and partial CSI, 
performance gain of the best optimized code over the standard 
code is about 0.4 dB, while for the no CSI case, improvement 
is about 0.3 dB. For M = 8 and partial CSI, improvement is 
0.4 dB, while no CSI shows 0.3 dB of improvement. 

V. Conclusion 

This work has presented an optimization procedure for 
FDPC Tanner Graph variable node degree distribution in 
the physical-layer network coded two way relay channel. 
The codes are optimized for relay reception in the multiple 
access phase. FDPC Codes from the DVB-S2 and WiMax 
standards are used to establish baseline error-rate performance. 
Optimized codes are generated under a variety of system 
configurations: modulation order M = 4 and M = 8, The 
optimized codes outperform the standard codes by margins 
from 0.3 to 1.2 dB. 
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